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ABSTRACT 


Deformation  under  conditions  of  uniaxial  strain  was 
found  elsewhere  to  be  very  nearly  recoverable  for  rocks  like 
granite,  diabase  and  certain  limestones.  This  apparent 
elastic  behavior  is  examined  here  in  greater  detail. 
Volumetric  compression  as  a  function  of  mean  stress  was 
closely  predictable  from  our  measurements  of  compressibility. 
However,  Poisson's  ratio  in  uniaxial  strain  was  up  to  25 
percent  higher  than  the  value  we  obtained  in  a  direct  static 
measurement.  Motion  on  cracks  was  assumed  to  have  occurred 
during  uniaxial  strain,  and  an  elastic  analysis  was  carried 
out  based  on  the  model  of  a  crack-filled  solid  developed 
previously  by  Walsh.  The  theory  correctly  predicted 
characteristics  of  volumetric  compression  and  Poisson's 
ratio  observed  experimentally. 

The  pressure  calculated  in  a  shock  experiment  for 
Westerly  granite  would  be  up  to  30  percent  higher  if  the 
appropriate  Poisson's  ratio  for  uniaxial  strain  is  used 
instead  of  the  commonly  tabulated  value  found  at  high 
hydrostatic  pressure. 
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INTRODUCTION 


Deformation  under  the  constraint  that  two  principal 
strains  remain  zero  is  called  uniaxial  strain.  Strain  is 
uniaxial  in  some  modes  of  plane  wave  and  shock  wave 
propagation  as  well  as  in  the  uniform  loading  of  bodies 
of  large  lateral  extent.  Both  situations  are  of  some 
geologic  interest;  shock  loading  of  rocks  is  used,  for 
example,  in  studies  of  the  equation  of  state.  It  is 
often  assumed  that  deformation  in  upper  parts  of  the 
earth  is  uniaxial. 

In  a  recent  laboratory  study  [1]  some  fifteen 
different  rocks  were  loaded  in  uniaxial  strain  to  obtain 
a  better  understanding  of  the  role  of  mineralogy  and 
porosity  in  their  deformation,  and  to  compare  any  failure 
with  that  in  conventional  triaxial  experiments.  The  rocks 
fell  into  two  groups,  those  of  low  porosity,  for  which 
deformation  even  to  very  high  stress  was  recoverable, 
and  those  of  high  porosity,  which  underwent  permanent 
compaction.  Behavior  of  the  second  group  is  described 
in  [1,  2],  Here,  we  explore  the  extent  to  which 
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the  recoverable  deformation  of  the  first  class  could  have  been 
predicted  from  independent  measurement  of  elastic  properties. 
Particular  attention  is  given  to  Poisson's  ratio,  for  this  is 
readily  obtainable  from  the  uniaxial  strain  experiment.  We 
use  here  the  'tangent*  value  of  Poisson's  ratio,  determined 
from  differential  changes  in  deformation,  rather  than  the 
'secant'  value.  Thus,  Poisson's  ratio  in  uniaxial  strain  is 
calculated  by  means  of  the  relationship  found  from  elastic 
theory,  v  =  da  3/ (da  1  +  da  3),  where  the  subscripts  1  and  3 
refer  to  the  axial  and  lateral  directions,  respectively. 
Similarly,  Poisson's  ratio  in  uniaxial  compression  is  -de 3/de 1 , 
where  de  is  the  differential  change  in  strain. 

Review  of  previous  work  and  many  experimental  details 
are  given  in  [1] .  The  method  used  in  [1]  was  nearly  identical 
with  that  of  Brown  and  Swanson  [3],  A  jacketed  cylindrical 
sample  was  so  loaded  by  axial  stress  and  confining  pressure 
that  the  circumferential  (and  therefore  the  radial)  strain 
measured  by  a  strain  gauge  remained  zero.  A  maximum  of  10  kb 
confining  pressure  could  be  applied;  at  this  limit  the  axial 
stress  reached  a  value  of  12  to  30  kb,  depending  on  rock  type. 
The  axial  stress  ox,  the  radial  stress  o3,  and  the  axial  strain 
£  1  were  recorded  during  the  experiment.  Permanent  strain 
following  a  cycle  of  loading  was  noted,  either  from  strain 
gauge  readings  or,  for  very  large  strains,  from  measurement  of 
the  sample  dimensions. 
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We  describe  here  additional  experiments  on  those  rocks 
which  showed  negligible  permanent  deformation.  These  tests 
were  carried  out  to  investigate  the  extent  to  which  recoverable 
behavior  indicated  true  elastic  deformation.  The  rocks,  total 
porosity,  and  modal  analysis  are  listed  in  Table  1.  Volume 
compressibility  and  Poisson's  ratio  were  determined  as  a 
function  of  confining  pressure,  using  the  same  samples  and 
strain  gauges  as  for  uniaxial  loading.  When  these  data  were 
compared  with  appropriate  parameters  from  uniaxial  loading, 
Poisson's  ratio  was  found  to  differ  significantly  from  that  of 
the  uniaxial  experiment.  We  suggest  below  an  explanation  for 
this,  based  on  microcracks,  following  analysis  developed  in 
earlier  studies  [4,  5]  of  the  effect  of  cracks  on  elasticity 
of  rock. 


2.  EXPERIMENTAL  OBSERVATIONS 

Volumetric  compressions  to  10  kb  pressure,  following 
the  procedure  outlined  in  [6],  are  listed  in  Table  2.  Volumetric 
compressions  have  been  assumed  to  be  three  times  the  measured 
axial  linear  compressions.  Where  the  rock  is  relatively  iso¬ 
tropic  this  will  be  close  to  the  actual  volumetric  compression; 
where  the  rock  is  anisotropic,  it  will  not  be,  but  it  is 
probably  the  most  appropriate  quantity  to  compare  with  ei  from 
the  uniaxial  strain  experiment. 
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TABLE  1.  ROCKS  STUDIED 


Rock 

Porosity 

% 

Modal  analysis 

Diabase,  II 

0.1 

49 

an45,  46  pyr,  3  ox,  2  mica 

Frederick,  Me.. 

Gabbro 

0.2 

70 

an^r  12  mica,  8  pyr,  7  am,  3  ox 

San  Marcos,  Cal. 

Schist 

0.3 

40 

qu,  26  mica,  15  or,  7  an5 ,  7  gar, 

Source  unknown 

5  ox 

White  marble 

0.3 

99 

ca 

Source  unknown 

Granite 

0.6 

26 

qu,  25  or,  37  anJ0,  12  mica 

Barre,  Vt. 

Granite 

0.9 

27 

.5  qu,  35.4  mi,  31.4  an17,  4.9  mica 

Westerly,  R.I. 

Abbreviations : 

qu 

quartz 

or 

orthoclase 

ca 

calcite 

pyr 

pyroxene 

an 

plagioclase 

mica 

mica,  clay 

gar 

garnet 

ox 

oxides 

mi 

microcline 
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TABLE  2.  VOLUME  COMPRESSIONS 
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Samples  of  five  of  the  rocks  were  set  up  just  as  for  the 
uniaxial  strain  experiment  II] .  An  axial  load  of  several 
kilobars  was  applied  at  a  number  of  different  confining 
pressures  starting  at  about  1  kb.  Both  axial  and  circumferential 
strains  were  observed,  enabling  ratio  to  be  determined  as  a 
function  of  pressure  to  10  kb.  The  values,  listed  in  Table  3, 
have  an  uncertainty  of  about  10  percent. 


3.  ANALYSIS 

Before  comparing  elastic  parameters  from  the  uniaxial 
strain  experiment  with  intrinsic  values,  we  consider  what  the 
former  ought  to  be ,  based  on  the  model  used  in  previous 
analyses  [5],  namely  a  homogeneous  elastic  matrix  containing 
an  isotropic  network  of  cracks.  We  will  consider  the  role 
these  cracks  play  in  material  loaded  in  such  a  way  that  strain 
is  uniaxial.  In  particular,  compressibility  and  Poisson's 
ratio  will  be  calculated  to  see  how  they  depend  on  crack 
parameters . 

The  effective  elastic  properties  of  the  model  under 
uniaxial  strain  are  found  conveniently  by  applications  of 
Betti's  reciprocal  theorem.  We  assume  only  that  the  matrix 
is  linearly  elastic  and  that  a  sufficient  number  of  cracks  of 
small  enough  size  exist  that  the  overall  deformation  appears 
to  be  homogeneous.  The  body  with  a  single  crack  or  section 
of  a  crack  acted  upon  by  the  axially  symmetric  stress 
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TABLE  3.  POISSON'S  RATIO 


distribution  (dot/  do3,  do3),  is  shown  in  Fig.  la.  To  find 
the  effective  lateral  strain,  we  consider  the  same  body  to  be 
loaded  on  both  external  and  internal  surfaces  by  the  stresses 
(0,  doU3»  doua),  as  shown  in  the  other  diagram  in  Fig.  lb. 

The  reciprocal  theorem  gives  the  equality 

doideUi  +  2do3deu3  =  2dou3de3  -  /  dxu3d6f  dA  (1) 

where  dxu3  is  the  shear  component  6f  dou3  in  the  plane  of  the 
crack,  and  the  integration  is  over  all  interior  surfaces.  We 
introduce  the  requirement  that  de 3  *  0  for  uniaxial  strain  and 
the  identities 


dTu3  =  sin6c°»0doU3 
<JeU3  -  U-u)  dou,/E 

deul  -  -2U  daUJ/E  (2) 

where  u  and  E  are  Poisson's  ratio  and  Young's  modulus  of  the 
matrix  material.  Equation  1  can  now  be  written 

do3/dOj  =  u/(l-u)  +  I  (3) 

where  2l(l-u)  -El  singcosg (d6^/do j ) dA 

Note  that  I  is  a  positive  quantity.  Brace  [1]  calculated 
Poisson's  ratio  from  the  slope  of  plots  of  a 3  versus  Ci  by 
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means  of  the  relation 


do j/do i  *  ugff/ (i  -  ufiff)  (4) 

where  is  the  effective  Poisson's  ratio  of  the  porous 

sample.  Combining  (3)  and  (4)  and  rearranging  gives 

°eff  "  0  +  <l-u) 1/(2  +  I)  (5) 

Thus,  u  f  must  be  greater  than  the  intrinsic  value  v  because 
the  second  term  on  the  right-hand  side  of  (5)  is  positive. 

The  stiffness  is  found  in  a  similar  way,  except  that  a 
uniform  axial  stress  doul  is  applied  as  in  Fig.  lc.  One  finds, 
following  the  steps  above,  that  the  effective  axial  compliance 
de,/do1  is 

del/do1  -  1/E  -  2udoJ/Edo1  +  I  (6) 

Introducing  (3)  into  (6)  gives 

de x /do i  -  (1  +  u)/3K(l  -  u)  +  I/3K  (7) 

The  first  term  on  the  right-hand  side  of  equation  (7)  is  the 
axial  compliance  of  a  non-porous  body  under  uniaxial  strain; 
the  second  term  I/3K  is  greater  than  zero,  so  the  effective 
compliance  is  increased  as  a  result  of  slip  on  interior  surfaces. 
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Note  that  dei  is  the  change  in  volumetric  strain  because  the 
other  components  of  strain  are  zero  here. 

An  expression  for  the  volumetric  strain  involving  the 
average  pressure  can  be  found  by  applying  the  reciprocal 
theorem  to  the  stress  state  in  Fig.  la  and  a  uniform  pressure 
dpu  applied  to  exterior  and  interior  surfaces  as  in  Fig.  Id. 
The  work  balance  from  the  reciprocal  theorem  is 

da i  deui  +  da3(deU2  +  deuj)  =  dp  (de^  +  2de 3)  (8) 

There  is  no  contribution  in  equation  (8)  from  slip  on  closed 
cracks  because  pressure  dpy  is  normal  to  this  displacement. 

If  cracks  are  distributed  isotropically, 

dEUl  ■  deU2  =  deUl  ‘  dP/3K  (9> 

Introducing  the  identities 


dV  =  dex  +  2de 3 
dpavg  =  (dai  +  2da3>/3 

we  find 

dV/Vdpavg  =  1/K  (10) 

Thus,  the  compressibility  determined  in  uniaxial  strain  when 
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all  cracks  are  closed  is  equal  to  the  intrinsic  value.  This 
result  can  be  verified  by  means  of  equations  (5)  and  (7) . 


4.  DISCUSSION 

The  theory  above  predicts  the  behavior  of  rocks  containing 
an  isotropic  distribution  of  cracks  under  uniaxial  strain. 
Poisson's  ratio,  when  applied  compressive  stresses  are  high 
enough  that  virtually  all  cracks  are  closed,  should,  based  on 
equation  (5) ,  be  greater  than  the  intrinsic  value  u  because  of 
the  relative  motion  between  crack  faces  which  may  occur  when 
applied  stresses  are  nonhydrostatic.  Poisson's  ratio  at  low 
stress  when  cracks  are  open  has  been  shown  theoretically  to  be 
less  than  the  intrinsic  value  u  [4],  Thus,  one  would  predict 
that  Poisson's  ratio  under  uniaxial  strain  for  crystalline  rock 
containing  cracks  should  increase  from  a  value  less  than  v  to 
a  value  greater  than  u  as  applied  stress  increases.  The 
magnitude  of  the  increase  cannot  be  evaluated  easily;  it 
depends  in  a  complicated  way  upon  the  number  of  cracks  and 
their  frictional  characteristics,  neither  of  which  are  now 
known  with  certainty  from  independent  measurements. 

We  test  the  above  prediction  for  Barre  granite,  diabase, 
gabbro,  and  marble  in  Figs.  2  and  3.  Apparently  the  values  of 
Poisson's  ratio  in  uniaxial  strain  (from  [1])  do  in  fact  exceed 
the  intrinsic  values  (from  Table  3)  for  marble,  granite,  and 
gabbro  at  high  stress.  At  low  stress  the  reverse  holds  for 
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gabbro,  granite,  and  diabase.  The  theory  cannot  be  tested 
against  the  data  for  diabase  at  high  stress  or  marble  at  low 
stress  because  the  differences  between  Poisson's  ratio  in  the 
two  tests  is  less  than  the  uncertainty  in  the  measurements. 

The  very  high  values  of  Poisson's  ratio  at  high  stress 
for  the  marble  in  uniaxial  strain  suggest  that  additional 
factors  than  motion  on  cracks  may  be  significant  for  this 
material.  The  most  likely  is  that  plastic  flow  of  calcite 
has  occurred;  this  would  probably  act  in  the  same  sense  as 
motion  on  cracks  and  cause  Poisson's  ratio  to  increase.  Some 
microscopic  evidence  presented  in  [1]  suggests  that  flow  of 
calcite  has  in  fact  occurred. 

One  can  carry  this  comparison  one  step  further  for 
Westerly  granite.  In  a  previous  analysis  [4]  for  the  condition 
of  uniaxial  stress,  it  was  also  found  that  Poisson's  ratio 
should  exceed  the  intrinsic  value;  comparison  of  that  analysis 
with  the  present  shows  that  Poisson's  ratio  under  uniaxial 
strain  should  be  less  than  the  value  under  uniaxial  stress, 
because  slip  on  cracks  is  inhibited  by  the  lateral  stresses 
in  the  former.  Thus,  Poisson's  ratio  for  the  proposed  model 
should  have  the  following  characteristics:  at  low  stress, 
Poisson's  ratio  is  less  than  the  intrinsic  value,  but  equal 
for  all  modes  of  deformation.  The  value  increases  as  stress 
increases  to  a  value  above  the  intrinsic  value.  Poisson's 
ratio  in  this  range  depends  upon  the  stress  state,  with  higher 
values  associated  with  conditions  of  less  lateral  constraint. 
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Measurements  on  Westerly  granite  seem  to  bear  out  tnis 
prediction.  Shown  in  Fig.  4  are  values  of  Poisson's  ratio  as 
a  function  of  axial  stress  under  conditions  of  uniaxial  stress, 
uniaxial  strain,  and  hydrostatic  pressure.  At  high  stress, 
Poisson's  ratio  in  uniaxial  strain  is  less  than  the  value  in 
uniaxial  stress  and  greater  than  the  value  under  hydrostatic 
pressure,  as  predicted  by  the  theory.  The  values  at  zero 
stress  are  less  than  the  intrinsic  value,  as  predicted,  although 
the  values  do  not  appear  to  be  equal.  This  discrepancy  may  be 
due  to  differences  between  individual  samples,  or  to  lack  of 
precision  in  the  measurements;  both  effects  are  more  significant 
at  low  stress  than  at  high  stress. 

Turning  next  to  compressibility,  we  find  that  the  theory 
above  predicts  that  volumetric  strain  per  unit  average  pressure 
is  equal  to  the  intrinsic  compressibility  for  any  state  of 
applied  stress  which  is  high  enough  to  close  all  cracks.  This 
equality,  which  is  a  property  of  an  ideally  elastic  material, 
was  unexpected  here  because  the  model  we  analyzed  is  not  truly 
elastic  in  this  region;  stiffness  in  the  axial  direction  is 
less  than  the  intrinsic  value  and  Poisson's  ratio  is  greater, 
for  instance. 

The  model  does  behave  elastically  at  low  stress  when  all 
cracks  are  open,  so  volumetric  strain  per  unit  average  pressure 
equals  compressibility  in  this  region,  as  well.  Compressibility 
is  less  than  the  intrinsic  value,  however,  because  of  the 
presence  of  open  cracks.  At  intermediate  stresses,  where 
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cracks  normal  to  the  maximum  compression  are  closed  and  those 
more  parallel  are  open,  the  rock  is  elastically  anisotropic, 
and  the  equality  does  not  hold.  Plots  of  average  pressure,  oj, 
versus  volumetric  strain,  ei,  under  uniaxial  strain  and  under 
hydrostatic  pressure  should  therefore  have  the  following 
characteristics.  The  curves  coincide  at  low  stresses  where 
all  cracks  are  open,  and  begin  to  diverge  as  stress  is  increased. 
At  high  stress  the  curves  are  parallel  but  not  necessarily 
coincident. 

In  Fig.  Z ,  measurements  from  [1]  are  compared  with 
intrinsic  volume  compressions  from  Table  2.  At  high  stress, 
data  in  uniaxial  strain  for  all  rocks  except  marble  lie 
approximately  parallel  to  the  corresponding  curve  from  tests 
under  hydrostatic  pressure.  Results  for  marble  deviate  farthest 
from  predicted  behavior,  perhaps  because  of  true  plastic  flow 
in  calcite  grains  which  may  have  occurred  in  this  experiment. 

In  previous  tests,  cracks  in  typical  rocks  seem  to  be  completely 
closed  by  pressures  of  2-4  kb,  and  most  of  the  effects  resulting 
from  their  closure  are  small  after  1  kb.  The  early  stages  of 
deformation  where  closing  of  cracks  is  important,  therefore, 
cannot  be  analyzed  in  plots  on  the  scale  of  that  in  Fig.  3. 

Note,  however,  that  the  separation  between  data  from  the  two 
tests  is  not  large,  suggesting  that  the  effect  of  anisotropy 
at  low  stress  upon  later  deformation  is  small. 

We  conclude  that  behavior  of  the  low  porosity  rocks  was 
not  strictly  elastic  in  the  uniaxial  strain  experiments  [1] , 
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even  though  they  did  appear  to  recover  in  an  elastic  fashion. 
Although  volumetric  compression  could  have  been  predicted  from 
independent  measurement  of  elastic  properties  (Fig.  5)  Poisson's 
ratio  could  not  (Figs.  2,  3,  and  4)  .  Presumably  otheir  elastic 
parameters  which  involve  deviatoric  stress  or  strain ,  such  as 
Young's  modulus  and  possibly  compressional  wave  velocity,  would 
have  shown  similar  disagreement. 

Nonuniqueness  of  Poisson's  ratio  has  important  bearing 
on  investigations  of  the  equation  of  state  of  rocks  because 
shock  data  is  obtained  from  experiments  which  apparently  are 
in  effect  uniaxial  strain  tests  at  high  strain  rate.  Poisson's 
ratio  must  be  known  in  order  to  reduce  this  data  to  obtain  plots 
of  average  pressure  versus  volume,  and  generally  a  value 
measured  in  tests  under  hydrostatic  pressure  has  been  used. 
Values  of  Poisson's  ratio  of  0.23  and  0.18,  corresponding  to 
measurements  at  high  and  low  pressure,  have  been  used  in 
reducing  data  on  Westerly  granite,  to  cite  only  one  example. 

The  average  pressure  is  15%  to  30%  higher  if  the  value  of  0.30 
obtained  in  uniaxial  strain  by  Brace  [1]  is  used. 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Body  containing  a  closed  crack  loaded  by  appropriate 
states  of  stress  for  applying  the  reciprocal  theorem. 

(a)  The  given  applied  stresses  (doi,  da3,  do3)  produce 
slip  of  d<5f  on  a  closed  internal  surface.  In  uniaxial 
strain,  d^a-O. 

(b)  A  uniform  stress  (0,  do3,  do3)  is  applied  to 
external  and  internal  surfaces  to  find  the  effective 
lateral  compliance. 

(c)  A  uniform  stress  (douj,0,0)  is  applied  to  external 
and  internal  surfaces  to  find  the  effective  axial 
compliance. 

(d)  A  uniform  pressure  dpu  is  applied  to  internal 
and  external  surfaces  to  find  the  volumetric  strain 
resulting  from  (dcfj.dc^,  da3). 

Poisson's  ratio  for  marble  and  Barre  granite  versus 
lateral  stress,  o3,  under  uniaxial  strain  (open  circles) 
and  hydrostatic  compression  (closed  circles) .  Error 
bars  show  the  typical  uncertainty  in  the  measurements. 
From  [1] . 

Poisson's  ratio  for  gabbro  and  diabase  versus 
lateral  stress,  o3,  under  uniaxial  strain  (open 
circles)  and  hydrostatic  pressure  (closed  circles) . 
Error  bars  show  the  typical  uncertainty  in  the 
measurements.  From  [1]. 
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Figure  4 


Figure  5 


Effective  Poisson's  ratio  o  for  Westerly  granite 
versus  axial  stress  cx  under  uniaxial  stress  [4], 
uniaxial  strain  [1],  and  hydrostatic  pressure  [1]. 

The  sample  used  in  the  uniaxial  strain  and  hydro¬ 
static  pressure  tests  was  different  from  the  sample 
used  in  the  uniaxial  stress  test.  The  error  bars 
show  the  typical  uncertainty  in  the  measurements. 

Volumetric  compression  compared  for  hydrostatic 
pressure  and  uniaxial  strain.  Curves  show  pressure, 

P,  versus  volumetric  strain,  0,  from  Table  2.  Circles 
give  mean  stress,  o,  versus  volumetric  strain,  ei, 
from  tests  under  uniaxial  strain  [1] .  Probable  error 
in  both  sets  of  measurements  was  less  than  the 
diameter  of  the  circles. 
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(c)  (d) 


Fig.  1  Body  containing  a  closed  crack  loaded  by  appropriate 
states  of  stress  for  applying  the  reciprocal  thcoror 
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Fig.  2  Poisson's  ratio  for  marble  and  Barre  granite  versus  lateral 
under  uniaxial  strain  (open  circles)  and  hydrostatic  compre: 
(closed  circles) 


Poisson* 


